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The interaction between static quarks is derived by applying many-body techniques 
i to QCD in Coulomb gauge. The result is shown to be exact in the IR and UV 

limits, and agrees remarkably well with lattice computations. 
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1 Introduction 



' One of the key impediments to understanding QCD at low energy is a quan- 

0^ . titative description of confinement. Here, we demonstrate how confinement 
may be constructed with a nonperturbative many-body analysis of QCD in 

. Coulomb gauged. 

' The starting point is the QCD JHamiltonian in Coulomb gauge, as derived 

, by Schwingercl and Christ and LeeoQ: 
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H = J dxtpi (-ia ■ V + 0m) ip + i J dx (^J- 1/2 UJ ■ UJ 

- g J dbap^a. ■ Atp + H c . (1) 
where Hq is the nonabelian version of the instantaneous Coulomb interaction: 

Hc = \J dxdyj-^p'ix) J 1/2 K ab (x, y; A)J 1/2 p b (y)J~ 1/2 . (2) 

In this expression the total charge density is denoted by p, J is the Faddeev- 
Popov determinant, J = dct(V ■ D), D is the covariant derivative in adjoint 
representation, D ab = <5 afc V — gf abc A c , and K is the nonabelian kernel 

K ab (x,y;A) = < x , a |^(-V 2 )^^|y, 6), (3) 

We shall make use of the inverse Faddeev-Popov operator, d = <?/(V • D) in 
the following. 
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Coulomb gauge is particularly efficacious in describing QCD in the con- 
finement region because all degrees of freedom are physical (no gauge con- 
straints need be applied) and the instantaneous interaction provides an ex- 
cellent starting point for the bound state problem. Furthermore, it is easy 
to see that the Coulomb interaction must yield confinement for static quarks 
because transverse gluon exchange is suppressed in the heavy quark limit. 
Thus it is clear that determining the behaviour of K is of central importance 
to developing an understanding of low energy QCD. 



2 IR and UV Behaviour of QCD in Coulomb gauge 

The evaluation of the Coulomb interaction is a nontrivial endeavour! However, 
the problem simplifies dramatically if one is willing to focus on the leading 
infrared (IR) behaviour of the theory (which, of course, is where confinement 
lies). To understand this statement it is useful to recall an old story in many- 
body physics. 

A naive attempt to compute the ground state energy of an electron gas 
(with an assumed positive background) at higher order is confounded by IR di- 
vergences arising from the instantaneous (abelian) Coulomb interaction. This 
problem was resolved many years ago by Bohm and PinesEl who showed that 
summing the leading IR-divergent diagrams (the ring diagrams) effectively re- 
places the divergent Coulombic 1/q 2 behaviour with a convergent l/q 2 e(to, k) 
behaviour. In fact, the sum has a wholly different analytic structure than any 
finite order in perturbation theory can produce. We now demonstrate that a 
similar situation occurs in the nonabelian case. 

The vertices of QCD in Coulomb gauge are shown in Figure 1 (there are 
additional vertices where n transverse gluons are emitted from the Coulomb 
interaction). We denote these vertices by a and call the number of ver- 
tices of a given type in a time-ordered diagram, n a . We anticipate that 
the Coulomb kernel K does indeed cause an IR divergence when inserted in 
a typical diagram. We denote the divergence due to if by A (for concrete- 
ness, one may consider A as an IR regulator). A general diagram then has 
an IR divergence of order A M where M = Dk — D g — D q . Here Dk, D g , 
D q , and D d represent the number of Coulomb, gluon, quark, and Faddeev- 
Popov propagators respectively. The following relationships may be easily 
verified!! 2D K = n ggK + n qgK + n Kg d, 2D d = n Kg d + 2n dgd + n qqd , and 
n-Kgd = n ggd + n qqd \ employing these leads to M = 1 + 2D K + D d - J2 a n a, 
and finally 
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M = 1 - n 4ff - n 3g - n qqg (4) 

This equation is of central importance - no diagrams are more diver- 
gent than K itself, and increasing the number of four gluon, three gluon, or 
quark-gluon vertices decreases the degree of divergence. Thus the leading IR 
behaviour of QCD is contained in matrix elements of the Hamiltonian itself - 
specifically it is not in iterations of H. 




Figure 1. Vertices in Coulomb gauge QCD. The solid line is a K propagator while the 
dashed line represents a d propagator. 



Before applying these results we considex the ultraviolet (UV) behaviour 
of QCD in Coulomb gauge. It is well knownEl that the beta function of QCD 
may be obtained directly from the Coulomb interaction (there is no need to 
consider loop corrections to vertices). Thus an accurate evaluation of K is 
sufficient to determine both the long distance and short distance behaviour of 
the static quark confinement potential! 

3 Coupled Gap Equations 

The previous observations simplify the evaluation of the Coulomb kernel: the 
leading IR divergent diagrams coincide with those of the rainbow ladder ap- 
proximation. Furthermore, vertex corrections to the rainbow ladder Dyson 
equation are suppressed in the IR limit (see Ref. [1] for more details). Thus 
K may be evaluated by solving a relatively simple nonlinear integral equation. 

An important additional step needs to be taken before evaluating K. 
Diagrams are all evaluated with respect to an implicit vacuum, normally taken 
to be trivial. However, a trivial vacuum is not likely to be a good starting 
point for describing confinement. Thus we choose to introduce extra freedom 
into the problem with the aid of a variational vacuum ansatz: 
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*o[A] = (A\uj) = exp 



~ J dkA a (k)cj(k)A a (-k) 



(5) 



Here, to is an unknown function to be determined by minimizing the vacuum 
energy density. 

The particular form of the trial vacuum is unimportant - it merely serves 
as a method to define the basis with which we choose to describe low energy 
QCD. The Fock space which is built on our variational vacuum consists of 
quasiparticles - constituent quarks and gluons. These degrees of freedom obey 
dispersion relations with infrared divergences due to the long-range instanta- 
neous Coulomb interaction of the bare partons with the mean field vacuum. 
This interaction makes colored objects infinitely heavy thereby removing them 
from the physical spectrum. However, color neutral states remain physical 
because the infrared singularities responsible for the large self-energies are 
canceled by infrared divergences responsible for the long-range forces between 
the constituents. This is the origin of the "colour confinement" phenomenon. 

Having specified the vacuum (up to the function w) we are ready to eval- 
uate the Coulomb interaction. Since the result depends on w we must also fix 
lo by minimizing the vacuum energy density (this is just the gap equation). 
However, the vacuum energy depends on {^q\Hc\^q), which depends on K\ 
Thus we must solve coupled equations for the gap and for K. Doing so yields a 
result for u> which may then be used to evaluate K - the equations themselves 
determine the gluonic (and hence quark) quasiparticle interaction. 




Figure 2. Comparison of an analytical approximation to ui(k) (dashed line) and the full 
numerical solution (solid line). 
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The details of this procedure are described in Ref. [1] . The result for lu is 
shown in Figure 2. The data are presented with respect to a scale m g = ui(0) 
which must be determined by comparison with experiment (recall that this is 
full QCD, the theory needs to be regulated and renormalizcd; the existence 
of m g is due to dimensional transmutation) . Finally, the Coulomb interaction 
may be evaluated. The result (after numerical Fourier transforming) is shown 
in Figure 3. Taking lattice data as "experiment" allows us to determine m g 
in terms of the lattice scale (the Sommer parameter). One finds m g ss 600 
MeV. 
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Figure 3. Static QQ ground state potential. The solid line is the numerical solution, data 
are taken from Ref. [8]. 

Finally, an examination of the UV limit of the coupled gap equations 
reveals that 12/11 of the one loop beta function is recovered. The reason 
is that the vacuum expectation value of the Coulomb interaction docs not 
include the gluon loop contribution to the running coupling. This is generated 
by iterating the Hamiltonian and may be incorporated if desired. Figure 3 
makes it clear that this will have little phenomenological effect. 

4 Conclusions 

Coulomb gauge allows dramatic simplification of the confinement (and bound 
state) problem of QCD. Confinement must be contained in the nonabelian 
Coulomb interaction. We have shown that the combination of allowing modest 
flexibility in the vacuum and a judicious evaluation of the Coulomb kernel is 
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sufficient to derive the lattice Wilson potential to great accuracy. 

Unfortunately, it is difficult to make cartoons which explain 'why' confine- 
ment occurs in this formalism. Nevertheless it is tempting to speculate that 
the origin of confinement is related to the Gribov problem - the nonuniqueness 
of Coulomb gauge. A simple practical resolution of the Gribov problem is to 
restrict the set of gauge configurations to those which form an absolute mini- 
mum of a functional in gauge transformation space (it is F[u] = Tr J d 3 x(A u ) 2 
where A" = uAw'-uVu^L The resulting set is called the fundamental modu- 
lar region. It may be showru that this region is convex, includes the origin, and 
intersects the Gribov horizon (defined as the location in gauge space where the 
Faddeev-Popov operator vanishes) . Two factors of the inverse Faddeev-Popov 
operator appear in the Coulomb interaction, and hence gauge configurations 
near the horizon are naively strongly suppressed. However the horizon forms 
a boundary in an infinite dimensional space and therefore carries a vast en- 
tropy. It is possible that these effects cancel, leaving a strong IR enhancement 
in the inverse Faddeev-Popov operator, d, and therefore in K. Details of this 
picture and phenomenological implications are being examined. 

Acknowledgments 

ES would like to thank the organizers of the Lake Louise Winter Institute 
for providing a truly marvelous venue for discussing physics and Dean Karlen 
and Francois Corriveau for organizing the first annual LLWI hockey game. 
This work was supported by DOE under contracts DE-FG02-00ER41135, DE- 
AC05-84ER40150 (ES), and DE-FG02-87ER40365 (AS). 

References 

1. A. P. Szczepaniak and E. S. Swanson, Phys. Rev. D 65, 025012 (2002). 

2. J. Schwinger, Phys. Rev. 127, 324 (1962). 

3. N. H. Christ and T. D. Lee, Phys. Rev. D 22, 939 (1980). 

4. D. Zwanziger, Nucl. Phys. B 485, 185 (1997). 

5. D. Bohm and D. Pines, Phys. Rev. 42, 609 (1953). 

6. A. R. Swift, Phys. Rev. D 38, 668 (1988). 

7. F. L. Feinberg, Phys. Rev. D 17, 2659 (1978). See also T.D. Lee, Particle 
Physics and Introduction to Field Theory, (Harwood Academic, New 
York, 1981). 

8. K.J. Juge, J. Kuti, and C.J. Morningstar, Nucl. Phys. Proc. Suppl. 63, 
326 (1998). 



llwi: submitted to World Scientific on February 1, 2008 



